Reliable Entanglement Detection Under Coarse— Grained Measurements 
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We derive reliable entanglement witnesses for coarse-grained measurements on continuous variable 
systems. These witnesses never return a "false positive" for identification of entanglement. We show 
that, even in the case of Gaussian states, entanglement witnesses based on the Shannon entropy can 
outperform those based on variances. We apply our results to experimental identification of spatial 
entanglement of photon pairs. 
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Introduction. The detection of quantum entangle- 
ment is crucial for the implementation of quantum in- 
formation tasks and technology. Typically, entanglement 
is detected through quantum state tomography or the 
measurement of entanglement witnesses, which involve 
fewer measurements [iHj]. For continuous variable sys- 
tems, quantum state tomography is difficult due to the 
large number of measurements required by the infinite di- 
mensional Hilbert space. Thus, entanglement witnesses 
are usually more appealing from an experimental point 
of view. However, even the experimental measurement 
of an entanglement witness for high-dimensional systems 
can be very time demanding. This could be due to the 
large number of measurements required to reconstruct 
the probability distributions associated to the continu- 
ous variables, for example see [Hj. It would thus be ad- 
vantageous to develop entanglement witnesses for coarse- 
grained measurements. This would allow for faster accu- 
mulation of experimental data with acceptable statistics, 
as well as for the use of less precise detection schemes. 

The usual entanglement witnesses for continuous vari- 
ables can fail under general coarse-grained measure- 
ments. In particular, improper application of a CV en- 
tanglement witness can result in a "false positive" in the 
case of extremely coarse-grained measurements. In this 
paper, using recently derived uncertainty relations 0,01, 
we develop new sets of entanglement criteria which are 
acceptable for measurements with any coarse graining. 
We illustrate the utility of our approach for spatial vari- 
ables of entangled photons 

Let us consider the Hilbert space H — Hi ® 712 of 
a continuous bipartite state pri- We have coordinate 
{xi,X2) and momentum {pi,P2) variables associated with 
the corresponding Hilbert spaces. Operators connected 
with these variables satisfy usual commutation relations: 
[xfe, Pj] = ihSkj, where k,j = 1, 2. 

We will consider entanglement witiresses involving the 
global variables: 



which obey the commutation relations: [x-i-,p_j_] = 2ih 
and [x±,p^] = 0. The marginal probability distributions 
for the bipartite state pi2 in the global variables picture 
(H]) read: 



R±ix±) = (a;±| / dx^ {x^\pi2\x^) ] \x±) , (2a) 



^± ip±) - {P±\ [^J dp^ {p^\pi2 Ipt)) \p±) , (2b) 

where \x±) and \p±) are eigenvectors of the operators x± 
and p_|_ respectively. 

A useful entanglement witnesses for continuous 
variables are the Mancini-Giovannetti-Vitali-Tombesi 
(MGVT) criteria which state that if pi2 is separable 
then (from now on we set fi. = 1): 

a^[R±]a^[S^] > 1. (3a) 

More sensitive tools are the entropic criteria Q which in 
the same situation provide the inequality 



h[R±] + h[S^] > ln(27re). 



X± = Xi ± X2, 



P± = Pi ± P2) 



(1) 



(3b) 



that is always stronger than ([5a|) except for the case 
of Gaussian states, when both criteria are equivalent. 
Here the variance [f] and continuous Shannon entropy 
h[f] of a probability distribution /(•) are defined in the 
usual manner: fT^[/] = / dz f [z) — (J dzzf{z)) and 
h[f] = — J dz f {z)ln[f {z)]. Since all separable states 
must satisfy inequalities pap and (|3b[) . their experimen- 
tal violation is an indication of quantum entanglement. 

Two couples: {R^,S+) and {R^,S^) of marginal 
probability distributions must satisfy the Hcisenberg 
uncertainty relation (HUR) and the Bialynicki-Birula- 
Mycielski (BBM) uncertainty relation however, with 
the same lower bounds as those in inequalities pap and 
((3b)) . Note that these lower bounds differ from the origi- 
nal HUR and BBM by the factor of 4 = 2^ and the term 
In 2 respectively. This is a consequence of Eq. (jlj , which 
describes the transformation {x+,x^) = U {xi,X2), with 
det U = 2. 
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A simple proof of witnesses ([3a|) and ()3bp can be based 
on the PPT criterion for positive partial transpose [lll - 
13 1 . If pi2 is separable then the operator p^^ (the oper- 
ator pi2 transposed with respect to the second Hilbcrt 
space %2) is positive. Thus, describes a physical 
state so that the HUR and BBM must be satisfied by 

^-R±^, j , as they are satisfied by {R±,S±). The par- 
tial transposition acts as a mirror reflection in phase 
space [isl , what implies invariance of the marginal posi- 
tion distributions, i.e. = R±, while the momentum 
distributions switch their roles since Sl^ = S^. This 
observation finalizes the derivation of both separability 
criteria. 

Entanglement witnesses under coarse graining. In 
practice, an experiment is performed with finite preci- 
sion. In the case of position and momentum of photons 
or massive particles this is due to the size of the detector. 
The experimental results are then discrete random vari- 
ables which represent the detection probability for each 
detection position. Let us assume that the continuous 
distributions are sampled at discrete points and let us 
introduce the detector function D, which we take to be 
the rectangle function: 

D^(z) = l^ forze [(j-i)?/,(j + i)ry] ^ 
1 elsewhere ' 

and define Zj = jrj so that D^{z)/ri — > (5(z — Zj) as 77 — > 0. 
In position and momentum spaces one shall obtain the 
detection probabilities {r^} and {s±} with values: 

00 CJO 

{r^}, = jdzDt{z)R±{z), {4}^ = jdzDf{z)S±{z), 

(5) 

related to fcth position and Zth momentum detectors re- 
spectively • The variances of these discrete prob- 
ability distributions can be obtained with 



k \^k) 



(6) 



and corresponding definitions for momentum variances 
CT^a . Points = fcA represent centers of the fcth detec- 
tor devoted to measurement of the x± variables. 

Discrete variances are good approximations to the vari- 
ances of the continuous distributions (f2a| and (j2b[) if the 
size A ((5) of the detectors is small relative to the prob- 
ability distributions R± {S±) Q. As A and 5 get large, 
the variances ct^a and a'^s tend to zero, since most of the 

continuous distributions R± and S± will eventually be 
localized in a single bin. However, the experimentalist 
is limited by the experimental precision of his/her de- 
tectors, which are on the order of A and 5. Thus, the 
variances calculated according to Eq. © are not valid 



estimators of uncertainty in the general case When 
the detectors happen to be too large, this can lead to a 
false detection of entanglement. For example, consider a 
pure Gaussian state, with momentum space wavcfunction 



'^{Pi,P2) = Aexp 



exp 



-{Pi-V2f 



(7) 

where A is a normalization constant. For (t+ =0"-, the 
state is separable. Nonetheless, if the size of the bins is 
about 3cr±, the MGVT criteria (|3a)) can be falsely vio- 
lated when the variance is calculated according to Eq. 

Instead of using the discrete probabilities directly, we 
can use them to construct approximations to the actual 
continuous probability distributions for the continuous 
variables x± and p±. Let us define the distributions: 



R^ (x±) 



00 



A;— — cx 



AT ^k 



sl{P±)= E {4} 

I— —00 



A 



DfiP±) 



(8a) 



(8b) 



so that R^ and Sj. go to R± and S± in the limit A,S 0. 
The continuous distributions R^ and arc the dis- 
cretized approximations to R± and S± obtained through 
coarse-grained measurements. Figures [2] and |3] show ex- 
amples of these continuous histogram functions. 

Calculating the variances of the distributions (|8ap and 
one has 0: [R^] = and [S^] = 

a'^s +6"^ /12, thus, the discrete variances in general under- 

estimate the inferred variances. As A and S grow large, 
these variances are given by A^/12 and (5^/12, which rep- 
resent the variances of the detector functions. 

In a similar fashion, the continuous Shannon entropies 
of the distributions (l8al) and (I8bl) are Idl: 



h[R^] = H [{r^}] + log A, h[Si] = H [{si}] + log 5, 

(9) 

where H [{r±}] and H [{s±}] are Shannon entropies 
corresponding to discretizations of the continuous distri- 
butions R± and S± [l3|, defined by the well known for- 
mula i?[{(7}] = — J2k Ik Ingfc. In the limit of large A, S 
the entropies ([9]) are given by the continuous entropy of 
the detector function (log A or log (5). 

Coarse-grained entanglement criteria. Generalizing 
the uncertainty relation established by Eq. (16) from 
Ref. (as in the previous case we shall multiply the 
lower bound by 4), wc obtain the inequality 



[R^] [Si] > 1 



(10) 



It has been shown (see an example given by Eqs. (17-19) 
from [^) that this relation is not simply a consequence of 
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FIG. 1: (color online) Experimental setup (left) and joint co- 
incidence distributions (right). Near- and far-field were mapped 
onto the detection planes by means of switchable lens systems (see 
text). The distance from the BBO crystal to the detection planes 
is 500mm. 



the HUR. However, ([TU]) becomes trivially satisfied when 
AS > 12. 

Moreover, we generalize the entropic uncertainty rela- 
tion given by Eq. (40) from Ref. (according to ([T|) we 
shall rescale A H' A/-\/2 and S H> S/y/^): 

^[{^±}] +^[{4}] >-ln(A,5)-ln[C(<5A)], (11a) 
with 



C(^)^nrin^^;i- 



(lib) 



Rao (^tV) UM denotes one of the radial prolate spheroidal 
wave functions of the first kind fisf . The lower bound in 
(lllap is nonnegative for all values of A, S. 

bmce and arc linear functionals of R± and S± 
we immediately generalize the previous PPT reasoning, 
i.e. if pi2 is separable then: 



a' [Rt] {S%\ - 1 > 



h\R^\ 



In [C ((5A)] > 0. 



(12a) 



(12b) 



Eq. p2ap is an entanglement witness based on the vari- 
ance product, while Eq. (|12bp establishes the entropic 
entanglement witness, both for the coarse-grained prob- 
ability distributions. In the limit A — >■ 0, (5 — >■ (|12ap 
goes to the MGVT criteria ([5a|) and (jl2bp reproduces 
the entropic criteria (j3b[) . Inequality (jl2bp is always 
stronger than (|12ap since even for a Gaussian quantum 
state, the coarse-grained probability distributions, illus- 
trated in Figures [2] and [3l for example, are not Gaussian. 

Experiment. We tested the coarse grained entangle- 
ment criteria using spatially entangled photons from 
spontaneous parametric down-conversion (SPDC), pre- 
pared approximately in the state ([7]). A 5mm long BBO 
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FIG. 2: (color online) Binned histogram distributions for ex- 
perimental results for a;_ measurements for different bin sizes 
n = 1,3,7,11. 



crystal was pumped with a 325nm cw pump laser in a 
TEMOO mode. Down-converted photons at the degen- 
erate wavelength of 650nm were detected through lOnm 
FWHM interference filters using single photon detectors. 
For a Gaussian pump beam, the transverse spatial struc- 
ture of the down-converted photon pairs can be approxi- 
mately described by a Gaussian wavefunction [1, lYil - 22 1 , 



which is factorable in the the x and y cartesian coordi- 
nates. Within this configuration, we use narrow slits in 
the detection system to access one of the transverse di- 
mensions of the two-photon field, whose continuous joint 
detection probability can be estimated from Eq. ([7]). 

As in Ref. [l^ , the spatial correlations were measured 
by using optical lens systems to map the near-field (x) 
and far-field (p) transverse coordinates onto the detec- 
tion planes. The x measurements used an imaging sys- 
tem with magnification of 4, consisting of a telescope 
with lenses fi ~ 50mm and /2 = 200mm, while the p 
measurements used a Fourier transform system with a 
lens /s — 250mm, as illustrated in Figure [TJ We mea- 
sured two-dimensional arrays of coincidence counts for 
the near-field (x) and far-field (p) variables by scanning 
the detectors across the vertical direction in the detec- 
tion planes. The widths of the slits used in our mea- 
surements were Sx = 0.050mm for x measurements and 
Sp = 0.020mm for p measurements. In both cases, the 
step size used in the scanning was equal to the slit width 
{sx or Sp). The measured joint probability distributions 
C(xi, X2) and C{pi,p2) are shown in the right-hand side 
of Fig. [1] These were then normalized to obtain proba- 
bility distributions i?(xi,a;2) and S{pi,p2). 

The marginal distributions i?^" and were calcu- 
lated from the joint distributions in Fig. [T] by group- 
ing the measurements into different size bins A„ = nAi, 
dm = rnSi, where n, m = 1, 3, 5, 7.... The smallest bin size 
is Ai = 2s^(/i//2) = 0.0250mm and (5i = 2sp(27r//3A) = 
1.546mm~^. Examples of the histogram distributions 
(before normalization) R^" and 5^'" are shown in Fig- 
ures [5] and 131 Figure U] (a) and (b) shows the entangle- 
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FIG. 3: (color online) Binned histogram distributions for ex- 
perimental results for measurements for different bin sizes 
m = 1,3,7,11. 
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ment witnesses (|12ap and (jl2b[) as a function of the num- 
ber of bin sizes n and to, respectively. The black region 
shows the area in which the witnesses do not detect en- 
tanglement. It can be seen that the entropic criteria iden- 
tifies entanglement for even larger bins than the variance 
criteria. This can be seen more distinctly in Figure |4] (c), 
which shows these same results for the case n = m. As 
can be seen, even for this approximately Gaussian state 
the entropic criteria outperform the generalized variance 
product criteria, due to the coarse graining. 

Error bars, which are smalUer than the symbols 
in Figure |4] (c), were calculated by propagating the 
Poissonian counts statistics, as well as the error in 
the center position of each bin, which we defined to 
be crA„ = 0.01%/2n/i//2 = 0.004mm and (7s„, = 
0.01^/2(2TO7r//3A) 0.55mm-i. Here 0.01mm is the 
minimum step size of the micrometers used to translate 
the detectors. 

Conclusions Building on previous results for entropic 
and variance-based uncertainty relations, we have pre- 
sented entanglement witnesses for continuous variable 
measurements that are obtained using discretized detec- 
tion systems. Our results have been tested in an exper- 
iment witnessing spatial entanglement in photon pairs 
obtained from parametric down-conversion. Due to the 
non-Gaussian nature of the binned histogram distribu- 
tions inferred from discretized measurements, the en- 
tropic entanglement witness performs better than the 
variance-based witness. These results should be appli- 
cable in other continuous variable quantum systems. 

Note added: Upon completion of this work, we be- 
came aware of similar results obtained by considering 
Einstein-Podolsky- Rosen-Steering inequalities [23 1. 
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FIG. 4: (color online) Evaluation of entanglement witnesses; (a) 
Entropic criteria l ll2bll and (b) Variance criteria Ill2all both as a 
function of the bin size n and m. Black denotes the region where 
criteria are not violated, (c) Evaluation of entanglement witnesses 
for the same bin size (n = m). The entropic criteria (blue squares) 
detects entanglement for a larger bin size than the variance criteria 
(red circles). 
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